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A STEFAN-TYPE STOCHASTIC MOVING BOUNDARY 

PROBLEM 


MARTIN KELLER-RESSEL AND MARVIN S. MULLER 


Abstract. Motivated by applications in economics and finance, in particular 
to the modeling of limit order books, we study a class of stochastic second- 
order PDEs with non-linear Stefan-type boundary interaction. To solve the 
equation we transform the problem from a moving boundary problem into a 
stochastic evolution equation with fixed boundary conditions. Using results 
from interpolation theory we obtain existence and uniqueness of local strong 
solutions, extending results of Kim, Zheng and Sowers. In addition, we formu¬ 
late conditions for existence of global solutions and provide a refined analysis 
of possible blow-up behavior in finite time. 
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1. Introduction 

Moving boundary problems allow to model multi-phase systems with separating 
boundaries evolving in time. Typically, the evolution of the free interface is strongly 
coupled with the evolution of the whole system. A classical example is the so 
called Stefan problem introduced in 1888 by Josef Stefan [Ste88], which describes 
the evolution of temperature v(t, x ) in a system of water and ice. In one space 
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dimension it reads as 


j%v(t,x) =rji-^ 2 v(t,x), x>x*(t), 

(!■!) %- t v{t,x) =r] w -j^ I v(t,x), x<x*(t), 

v(t,x*(t)) = 0, 

where r] w and rji are the thermal diffusivities of ice and water, and a;*(f) is the 
position of the interface between the two phases. The evolution of the interface is 
governed by the so-called Stefan condition 

(i.2) = e -- ^w(M*(*)-)), q > o. 

This problem and various extensions have been studied extensively in the second 
half of the 20th century, see [Vui93] for a review of the literature. For classical 
solutions of semi-linear extensions of (1.1) see e. g. [FP79], [Lun04]. In addition to 
the theory of classical and weak solutions, the corresponding evolution equations 
have been studied in the framework of maximal L p -regularity, see [EPS03], [PSS07] 
and references therein. Compared to the deterministic case, stochastic partial dif¬ 
ferential equations with free or moving interface have received much less attention. 
One exception is [BDP02] where Barbu and da Prato show existence of a solution 
and an invariant ergodic measure for the linear problem (1.1) with additive noise 
in multiple dimensions. 

More recently, both deterministic and stochastic moving boundary problems have 
been applied in economics and finance to dynamic models of trading, in particular 
to models of so-called (electronic) limit order books where orders of buyers and 
sellers participating in stock exchanges are stored, see e.g. [LL07, Zhel2, MTB14, 
BHQ14]. In such models, the space coordinate x typically corresponds to price 
(usually on logarithmic scale), and the quantity v(t,x) to the density of buyers or 
sellers willing to commit to a transaction at time t for the price x. Buyers are 
recorded with positive sign and sellers with negative sign, such that the two phases 
of the system distinguish buyers from sellers. Of particular interest is the evolution 
of the separating boundary, which corresponds to the marginal price at which both 
sellers are currently willing to sell and buyers are willing to buy. Zheng [Zhel2] for 
example proposes the following stochastic moving boundary problem as a model 
for dynamic trading in a limit order book: 

§i v {t,x) = r)s-£iv(t, x) + cr s (|a; - x*(t)\)d£ t (x), 

( L3 ) = Vb-§^v(t,x) + <r b (\x - x*(t)\)d£ t (x), 

v(t,x*(t)) = 0, 

where the subscripts b and s correspond to buyer and seller respectively and d£t {x) is 
Gaussian noise. The evolution of the interface is governed by the linear Stefan condi¬ 
tion (1.2). The accompanying mathematical theory is developed in [KZS12, KMS10] 
and numerical analysis in [KS12]. Other examples can be found in Lasry and Li¬ 
ons [LL07] where a free boundary model for price formation under negotiation is 
introduced in a mean-field game setting. Another part of the literature derives 
SPDE models for limit order books as functional limits from discrete queuing mod¬ 
els of orders that arrive and then are filled or cancelled. For examples of this 
approach see e. g. [BHQ14] where a parabolic SPDE as a model for the order book 
is obtained in the limit. In addition, there is a series of papers by Bouchaud et 
al. [MTB14], [DBMB14] with PDE and SPDE models observed as limiting equa¬ 
tions of particle models. 


x > x*(t), 
x < x*(t), 
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In this paper we study a Stefan-type stochastic moving boundary problem, which 
can be considered an extension of (1.3) and of the theory developed in [KZS12] with 
several important differences in scope and methodology: 

> Instead of the homogeneous linear stochastic Stefan problem, we allow for 
a more general drift coefficient and in particular a non-linear boundary 
condition replacing (1.2). Recent empirical studies of the dependency of 
price change on the imbalance of the order book (see [CKS14] and [LPS13]) 
suggest linear behaviour for balanced order books and non-linear behaviour 
when imbalance is large. 

> In addition to mild and weak solutions as in [KZS12] we obtain solutions 
in the analytically strong sense and make the transformation from free to 
fixed boundary, that is introduced in a deterministic setting in [Lun04] and 
used in [KZS12], rigorous in a stochastic setting. 

> We combine tools from the SPDE framework of da Prato and Zabczyk 
(cf. [DPZ14]) with results from interpolation theory, which allows for greater 
generality and avoids direct computations using the heat kernel as in [KZS12] 

2. A STOCHASTIC MOVING BOUNDARY PROBLEM 

2.1. Problem formulation. Our goal is to establish a framework for solving sto¬ 
chastic moving boundary problems of the type 


x > x*(f), 


with the moving boundary x* ( t ) governed by 

J = e( ■§£v(t,x*(t)+),£v(t,x*{t)-)'j, 
for t > 0, x G R, with Dirichlet boundary conditions at x*, i. e., 

v(t,x*(t)+) = 0, 
v(t,x*(t)~) = 0, 


(2.1) 









d v(t, x) = 

r l+ dx 2 

V + H+ 

(x- 

- X*(t),v, 

_d_ 

dx 

V ) 

d t 






+ <J+ (x 

- 

X* 

(t),v) 

d£t(x), 

d v(t, x) = 

n -£L 
dx 2 

V + fl- 

(x- 

-X*{t),v, 

d_ 

dx 

V \ 

d t 






+ (J- {x 

- 

X* 


d £t(x), 


(BC) 


X € 


for t > 0. The coefficients are functions /i± : R 3 —» R, a± : R 2 —>• R, and real 
numbers r]± > 0. We denote by £ the spatially colored noise given by 

(2.2) £ t (x) := [ T c dW s (x), T ( w{x) := [ C(x,y)w(y) dy, 

for some integral kernel £ : R 2 — > R and a cylindrical Wiener process W on the 
Hilbert space U = L 2 (R) with covariance operator identity. As usual, W lives 
on a filtered probability space (f2, jF, (JFt)t> o,P). For each t > 0 we require that 
x i-a v(t,x) is continuously differentiable on (—oo ,x*(t)) as well as on (x*(t),oo), 
such that all first derivatives appearing in (2.1) can be understood in the classical 
sense. The second derivative should be considered a weak derivative and a suitable 
function space for v as well as the precise notion of ‘solution’ to (2.1) will be defined 
below. 


We now make precise what we understand by a solution to the stochastic moving 
boundary problem (2.1). In general, solutions to the moving boundary problem 
may be local, i.e. only exist up to a stopping time r. To formalize this, it will be 
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convenient to work with stochastic intervals. Given two stopping times c < t the 
stochastic interval [?, t] is defined as 

[?, t] := {(u,t)e(lxR + : ?(w) < t < r(w)} . 

By using strict inequalities also the open stochastic interval ]?, t|[ and half-open 
analogues can be defined. As usual in a probabilistic setting we will soon ‘drop the 
omega’ and say e.g. that for two stochastic processes X and Y the equality X t = Y t 
holds for all t £ [0,r], when we mean that X t (co) = Y t (uj) for P-almost all ui and 
all t such that (w,f) £ [0,tJ, i. e. 

P [Xt = Yf,Vt < t\ = 1. 


To formalize the moving frame for the moving boundary problem we define for 
i£l the function space 


(2.3) T(x) := jw : R —>• R : v\^ {x} £ H 2 (R \ {x}) n H'(R \ {x})} , 


where Hq and H 2 are the usual Sobolev spaces. Note that due to the Sobolev 
embeddings, any function v in T(x) can be identified with an element of L 2 (R). 

Finally, using the notation from (2.1) we introduce the functions ~p : R 4 —> R, 
a : R 2 —> R, 


(2.4) 

Ji(x,v,v f ,v") 


V+v" +n+{x,v,v'), 
r\-v" + fi-(x, v, v'), 


x > 0, 
x < 0, 


cr(x, v ) 


cr + (x,v), 

a-(x,v), 


x > 0, 
x < 0. 


Definition 2.1. A local solution of the stochastic moving boundary problem (2.1) 
on the stochastic interval [0,t[[, with initial data vo and Xq, is a couple (u,x„) of 
stochastic processes, where 

(v,x*) : [0,t[—X U (r(x) x {x}) C L 2 (R) x R. 


such that (v,x*) is predictable as an L 2 (R) x R-valued process, and 

v{t)-v o= / p(. -x*(s),v(s), -j%zv(s)) ds 

Jo 

+ / a{. -x»(s),u(s)) d^ s (.), 

^x»(t) = e{^ i(t,x*(t)+), £v(t,x*(t )-)), 
x*(0) = xo 

holds on [0, rj. The first equality is an equality in L 2 (R); the first integral is a 
Bochner integral in T 2 (R), and the second one a stochastic integral in L 2 { R). 

The solution is called global, if r = oo and the interval JO, r[ is called maximal if 
there is no solution of (2.1) on a larger stochastic interval. 


2.2. Assumptions and main results. We introduce the following assumptions 
on the coefficients appearing in (2.1) 


Assumption 2.2. The functions are continuously differentiable and 
(i) there exist a £ L 2 (R+), b,b£ L“ C (R 2 ; R) such that for all x,y, z £ 


W±(x,y,z)\ + \ -§zH±{x,y,z) \ < a(|x|) + b{y , z) (\y\ + \z \), 


and 


-§^y±(x,y,z) + \£n±(x,y,z)\ < b(y,z), 


(ii) y± and their partial derivatives (in x, y and z) are locally Lipschitz with 
Lipschitz constants independent of x £ R. 
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Assumption 2.3. The functions a± are twice continuously differentiable and 

(i) For every multi-index I = (i,j) £ N 2 with |/| < 2 there exist aj £ L 2 (R+) 
and bj £ L“ C (R, R + ) such that 

ai{\x\) +6/(y) \y\ , j = 0, 
bi(y), j 7 ^ 0 . 


dx z dy° 


r(x,y) 


(ii) <r± and their partial derivatives (in x, y and z) are locally Lipschitz with 
Lipschitz constants independent of x £ R. 

(iii) <r± satisfy the boundary condition 

(2.5) cr±(0,0) = 0. 


Remark 2.4. Later on, certain Nemytskii operators will be defined through fi± and 
<j± and the assumptions made above can be traced back to requirements on the 
regularity of these operators, see Appendix A and B. Also note that if <r + or er_ 
is independent of x £ R, then, in Assumption 2.3, part (iii) is a consequence of 
part (i). 

Assumption 2.5. g : R 2 — > R is locally Lipschitz continuous. More precisely, for 
all N £ N there exists an L e ,iv such that 

I e{y) ~ Q(y)\ < L e , N I y - y\ for all \ y \, \ y \ < N. 


Assumption 2.6. £(.,y) £ C 3 (R) for all y £ R and ^y T £(x,.) £ L 2 (R) for all 
x £ R, i £ {0,1, 2, 3}. Moreover, 


( 2 . 6 ) 


sup 

aj£]R 




L 2 (R) 


< OO, 


= 0,1,...,3. 


For the rest of this paper, we use the notation := 

Example 2.7 (Convolution). Let C be a convolution kernel, i. e. C,(x 7 y) := — y ), 

x, y £ R. If C £ C'°°(R) Cl H 3 ( R), where H 3 denotes the Sobolev space of order 
3, then Assumption 2.6 is satisfied. In this case, the operator corresponds to 
spatial convolution with £. 


Example 2.8 (Stochastic Stefan Problem). Let = /r_ = 0, a+(x,v) = cr-(—x,v) = 
v and g(x±,X 2 ) = Q ■ ( x^ — xi) for some g £ R. Then, (2.1) is the two-phase Ste¬ 
fan problem with multiplicative colored noise. With r \2 = 0, = 0 and 

C{x,y) := £(x — z) we end up with the one-phase system discussed in [KZS12]. 
Even though our assumptions and proofs are formulated for the two-phase case, it 
is straight-forward to adapt them to a one-phase setting. 

Example 2.9 (Two-Phase Burger’s equation). The case 

fi+(x, v, v') := /j,-(—x, v , v') := v ■ v', x £ R>o, v, v' £ R, 

yields a stochastic version of a two-phase viscous Burger’s equation in one dimen¬ 
sion. Obviously, Assumption 2.2 on is satisfied. 

Example 2.10 (Reaction-Diffusion-type drift). Set p±(x,v,v') := f±(v), for some 
f± £ C 1 (R) with locally Lipschitz derivative and /+(0) = /_(0) = 0. Also in this 
case it is easy to check that Assumption 2.2 is satisfied. 


Let us also remark here that without substantial change in our proofs the con¬ 
stant Laplacian terms r]±-^jv in (2.1) can be replaced by space-dependent Lapla- 
cians in the divergence form (j]±(x — x*(f)) • -§^v) for some scalar functions r]± 
that are bounded by strictly positive constants. 


Our first main result concerns the existence of a maximal local solution to the 
moving boundary problem (2.1). 
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Theorem 2.11 (Maximal Local solution). Let Assumptions 2.2, 2.3, 2.5, and 2.6 
hold true and let xq G M and vq G T(a;o). Then there exists a predictable, strictly 
positive stopping time r and a local solution (v,x*) of (2.1) on the maximal interval 
JO, r[ in the sense of Definition 2.1. For almost every uj £ 12 it holds that v(lo; .) £ 
C([0, t(w)); H 1 (M.)) and x*( uj,.) £ C' 1 ([0, t(w)); K.). Moreover, (v,x*) is unique 
among all H 1 © R- continuous solutions. 

Remark 2.12. The continuity statement implies that for all x G R, t H > v(t,x) is 
almost surely. 

Imposing some additional assumptions on o± and p the solution becomes global. 

Assumption 2.13. The functions b and b in Assumption 2.2 are (globally) bounded, 
and there exist functions cr± G JL 2 (R+)nC 2 (R>o) and cr± G BUC 2 { R>o), the space 
of all functions with bounded uniformly continuous second derivative, such that 

<x+{x,y) = 0+0*0 + v 2 +(x)y, = P-(x) + a 2 _(x)y 

for all x G R>o and y G R. 

Theorem 2.14 (Global Solution). If p is bounded and Assumption 2.13 holds in 
addition to the assumptions of Theorem 2.11, then t = oo almost surely, i. e. (2.1) 
has a global solution. 

In Theorem 4.5 we provide a refined analysis of the case of finite-time blow-up 
(r < oo). It turns out that under Assumption 2.13, but with p unbounded a finite¬ 
time blow-up of the system (2.1) must coincide with a blow-up of the boundary 
terms. 


2.3. Overview of the proof. Our treatment of equation (2.1) consists of three 
steps 

> Transformation into an equation with fixed boundary; 

> Formulation of the transformed equation as an abstract stochastic evolution 
equation; 

> Solving the abstract evolution equation by a fixed-point-argument. 

For the first step we apply a change of coordinates 

(2.7) zii(f, x) := v(t,x*(t) + x), and u 2 (t, x) := v(t, x*(t) — x), 

i.e. new coordinates are defined relative to the free boundary x*(t), which yields 


du 1 (t,x) = + p + ( 2 c,ui, ^«i) 

d t + a+ ( x , ui) d£ t (x* ( t ) + x ), 




( 2 - 8 ) d u 2 (t,x)= r] 2 -^u 2 + p- (-x,u 2 ,-^u 2 ) 

d t + (J_ (—x,u 2 ) d£ t (x*(t) - a;), 


■(*<«) 


dX*(t) = e(-^u 1 (t,0+),-^u 2 (t,0+)) d t, 
for t > 0 and x > 0 with Dirichlet boundary conditions, 


(D) m(f, 0) = 0, u 2 (t, 0) = 0. 

Note that the classic chain rule is not sufficient to derive (2.8) from (2.1), since v is 
not differentiable in time. Rather a special case of Ito’s formula (a ‘stochastic chain 
rule’) is needed to justify the computation. The transformation turns the moving 
boundary into a fixed boundary at x = 0, but introduces an additional non-linear 
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and unbounded drift term involving the spatial derivatives -§^u\ and -§^u 2 . 

The second step is the abstract formulation of (2.8) in terms of the stochastic 
evolution equation 

(2.9) dX(t) = [AX(t) +B(X{t))} dt + C(X(t)) dW t , A(0) = X 0 . 
where 

X{t) := (m (i,. ),u 2 {t , ,),x*(t)), t > 0 

and W is a cylindrical Wiener process with covariance operator Id on the separable 
Hilbert space U = L 2 (R). Introducing the shorthand 

(2-10) i( M ) = (^ Wl ( 0+ ),^ U2 (0+)) 

for the boundary terms, the coefficients of (2.9) are given by 

(V+ A 0 0\ 

(2.11) A=\ 0 r/_ A 0 -cid, 

\ 0 0 0 / 

/ n + (x,U!(x), -§jUi{x)) + ^Ui{x) ■ g{l(u(t))) \ 

(2.12) B{u)(x) = \ n_(-x, u 2 (x), ^u 2 (x)) - ^u 2 (x) ■ + c id, 

V 6(Z(um ) 

( a + {x,ui(x))T^w(u 3 +x) \ 
a-(—x,u 2 (x))T^w(u 3 — x) ] , uG D(A), w £ U, x > 0. 

0 J 

Here, A is the Laplacian on R + with Dirichlet boundary conditions and c > 0 is 
an arbitrary constant, whose sole function is to move the spectrum of A into the 
negative half-line (—oo, 0). 


Finally, the solution of (2.1) in Theorem 2.11 will be obtained from the unique 
strong solution A' on [0,r[ of the stochastic evolution equation (2.9) with initial 
data X 0 := (n 0 (. + Zci)|R + ,i’o(a;o ~ Ok+^o) by setting X = {u 1 ,u 2 ,xA) and 


j(t,x) := 


Ui(t, x — £*(t)), x > x*(t), 
u 2 (t, x*(f) — x), x<x*(t), 


In the remainder of the paper we will make the above steps rigorous, by traversing 
them in the reverse direction: 

> In section 3 we show that under certain assumptions the abstract stochastic 
evolution equation (2.9) has a unique strong solution. 

> In section 4 we show that the parameter assumptions made in Section 2.2 
are sufficient for the assumptions of section 3 

> In section 5 we show the stochastic chain rule that is necessary to make the 
transformation to fixed boundary rigorous and collect all pieces to complete 
the proof of our main results. 

Remark 2.15. Considering carefully our proof of the existence result in the next 
section it can be seen that equation (2.9) can be solved also for homogeneous 
Neumann or even Robin boundary conditions. Of course the boundary conditions 
on a± in Assumption 2.3 (iii) have to be adapted accordingly. The main difference 
to the case of Dirichlet boundary conditions, is that a discontinuity at the boundary 
introduces a jump to the dynamics of v in equation (2.1) at any given point i£l, 
every time the boundary x*(t) crosses x. In particular the ‘stochastic chain rule’ 
developed in Section 5 is no longer sufficient to pass from the moving boundary 
equation (2.1) to the fixed boundary equation (2.8) and back. 
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Remark 2.16. The existence result for the centered equations (2.9) can be extended 
to the case with Brownian noise in the boundary, without any problems. That is, 

dx*(f) = D(u(t)) dt + cr* d B t 

where D is any locally Lipschitz operator from T>{A) into R, and <r* > 0. Here, B 
can be either independent of W, or a Hilbert-Schmidt transformation of W into R. 

3. Solving a stochastic evolution equation 

3.1. Preliminaries. In this section we concentrate on the evolution equation (2.9), 
i.e. 

(3.1) dX(t) = [AX(t) + B(X(t))]dt + C{X(t)) dW t , t> 0, 

where IT is a cylindrical Wiener process with covariance operator Id on a separable 
Hilbert space U . At this point it is sufficient to assume that X takes values in 
an arbitrary separable Hilbert space E with norm ||.||. On the coefficients A, B , C 
we will impose assumptions that are milder (but also more abstract) than the as¬ 
sumptions made in section 2 on the coefficients of the free boundary problem. As 
will be shown in section 4 the assumptions below are implied by the assumptions 
from section 2.2 such that eventually the results on the evolution equation (2.9) can 
be used to solve the free boundary problem (2.1). Nevertheless, the results of this 
section may be of independent interest when generalizations of (2.1) are considered. 

On the operator A in (3.1) we make the following assumption. 

Assumption 3.1. A is a densely defined and sectorial operator with domain 
T>(A) C E. Moreover, the resolvent set of A contains [0, oo) and there exists a 
M > 0 such that the resolvent R( A, A) satisfies 

(3.2) ||i?(A,A)|| < for all A >0. 

Remark 3.2. This assumption is equivalent to each of the following statements 

> Equation (3.2) holds and the resolvent set of A contains 0 and a sector 

{A £ C : | arg A < 9} 

for some 9 £ (7t/2,7t). 

> The operator A is sectorial and —A is positive in the sense of [Lun09|. 
Assumption (3.1) ensures that A generates an analytic semigroup (St)t> o and that 
suitable interpolation spaces between E and T>(A) can be defined through fractional 
powers of —A. They also imply that the semigroup St is of strictly negative type, 
i.e. there exist 5, M > 0 such that US'*|| < Me~ St . Note, that if M = 1 then St is 
a contraction semigroup, which we shall not assume a priori. 

Using the semigroup St that is generated by A we can introduce the important 
concept of mild solutions. 

Definition 3.3. Let X = (X(t)) be a 2?(A)-valued predictable process and t be a 
predictable stopping time. 

> X is called global mild solution to the stochastic evolution equation (3.1) 
on T>(A) with initial data Xq £ 2?(A), if 

X(t)=S t X 0 + f S t - s B(X(s)) d 3 + f S t - s C(X(s)) dW s . 

Jo Jo 

holds for all t > 0, P-a.s. 

> X is called mild solution on [0, t|[, if (3.1) holds on the stochastic interval 


(3.3) 
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> The stochastic interval [0,r[ is called maximal for X if there is no 'D(A)- 
continuous extension of X to a larger stochastic interval. 

In the last two terms of (3.3), f denotes the Bochner and stochastic integral on 
the Hilbert space T>(A), respectively. If we want to emphasize the underlying space 
T>(A) we write global mild l?(A)-solution and mild 2?(A)-solution respectively. 

Finally we will be able to show that the mild solution is also a strong one in the 
following sense: 

Definition 3.4. Given T>(„4)-valued initial data Xq and a predictable stopping time 
r, X is called strong solution of (2.9) on [0, t|, if X is a 2?(„4)-valued predictable 
process and 

AX(s) + B{X{s)) d s+ f C(X(s)) dW s , 

Jo 

holds on [0, tJ. Global solutions and maximality are defined in the same way as 
for mild solutions. 

3.2. Interpolation spaces. By taking fractional powers of —A we introduce inter- 
and extrapolation spaces for E. For a > 0 we define 

(3.5) E a := V((-A) a ), \\h\\ a := \\(-A) a h\\ E , h £ E a . 

It is known that also E a with the induced scalar product is a separable Hilbert 
space. In particular, ||.||^ is equivalent to the graph norm of A and the following 
continuous embedding relations hold for a £ [0,1]: 

(3.6) V(A) = E 1 ^E a ^>E 0 =E, 

Note that the restriction of A to any E a , a £ [0,1] is again a densely defined and 
closed operator on E a . Moreover, it is the infinitesimal generator of the restriction 
of St to E a , which is again an analytic (contraction) semigroup; see e.g. [EN00, 
Ch. II.5]. 

The following regularity property of St between different interpolation spaces 
E a , a £ [0,1] will be crucial in the proofs that follow. We derive it from results 
in [Lun09] on interpolation spaces. 

Lemma 3.5. Let /3 > 0 and a > (3. Then, for all t > 0 and h £ Ep, 

\\Sth\\ a <K a , p ^e- St \\h\\ fl . 

Note that the factor in front of ||h||^ is integrable at time t = 0, which is the key 
property used in the estimates concerning the mild formulation of (3.1) on E\. 

Proof. Suppose first that a = /? + n for some n £ N, then we get from [Paz92, 
Thm 1.5.2d and p. 70] that there exists K n > 0 such that 

(3.7) \\S t h\\ a = ||(- J 4)"S t (-^/»|| 0 < K n r n e~ st \\h\\ 0 . 

Now assume that a £ (/3 + n, f3 + n + 1) for some n £ No and set 0 = a ~ /? — n£ 
(0,1). By [Lun09, Prop. 4.7] the real interpolation space ( Eo , Ei)g ,i is continuously 
embedded into T>{{—A) 0 ). Combining this fact with [Lun09, Cor. 1.7] we obtain 
that there exists K > 0 such that 

\\h'\\ e < K H/i'll l~ 9 Wh'Wl for all h' € E x = T>(A). 

Now let h £ T>({—A)P) and set h' = (— A) n S t {—A)^h £ T>(A). Applying the above 
inequality and using boundedness of the semigroup St we obtain 


(3.4) X(t) =X 0 + f 

Jo 
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Finally, (3.7) for n and n + 1 yields 

\\S t h\\ a < K"e~ 5t t -^«-«("+ 1 ) \\h\lp , 

proving the result. □ 

To deal with the singularity in 0 on the right hand side above, we will use an ex¬ 
tended version of Gronwall’s lemma, see [Lun95, Lem 7.0.3] or, for a proof, [Hen81, 

p. 188], 

Lemma 3.6 (Extended Gronwall’s lemma). Let a > 0, a, b > 0, T > 0, and 
u : [OjT] — > R be non-negative and integrable. If, for all t £ [0,T], 

(3.8) u(t)<a + b ( u(s)(t — s)“ -1 ds, 

Jo 

then exists a constant K a ^^T, depending only on a, b and T, such that, 

(3.9) u(t) < aK at b'T, f£[0,T]. 

3.3. Existence of global mild solutions. We start by discussing global solu¬ 
tions. Subsequently, the existence of local solutions under milder assumptions will 
be shown by localizing with appropriate stopping times. Denoting by HS(t/, E\) 
the (Hilbert) space of Hilbert-Schmidt operators from U to E\ we introduce the 
following Lipschitz-type assumption, which will imply the existence of global mild 
solutions to (3.1). 

Assumption 3.7. There exists a £ (0,1] such that B : E\ — > E a and C : Ei —> 
HS(f7, E\) are Lipschitz continuous, i.e. there exists a constant L such that 

(3.10) \\B(Y)-B(Z)\\ a + \\C(Y)-C(Z)\\ 

HS(U,Ei) — LWY-ZW,, 

holds for all Y, Z £ E\. 

Remark 3.8. Assumption 3.7 implies a linear growth bound on B and C in the 
sense that 

(3.11) \\B(Y)\\ a + ||c(y)|| HS( ^ i) < m( i + Hry. 
for all Y £ Ei. 


Theorem 3.9 (Global Mild Solution of (3.1)). Let Assumption 3.1 and 3.7 hold 
true and let p > 1. Then, for every initial data Xq £ L 2p (fl;Ei) there exists a 
unique global mild solution X of (3.1) on E\. Moreover, 


(3.12) 


E 


sup \\X(t)\\l P 
0 <t<T 


<K (l + E 



for all T > 0 and X is E\-continuous almost surely. If, in addition, St is a 
contraction semigroup, then the statement is true even for p = 1. 


Remark 3.10. Without much effort one can extend the theorem to time- and path- 
dependent predictable coefficients B : Six R>o x E\ —> E a and C : fix R>o x£i -> 
HS([7, Ei), provided that (3.10) and (3.11) hold. The same is true for Theorem 3.17 
below. 


Proof. The theorem will be shown using a fixed-point argument. Using Lemma 3.5 
we will be able to prove that the following mapping is a contraction. We define 

K(Y){t) := S t X 0 + f S t - s B(Y(s)) ds + [ S t - B C(Y(s)) dW s , 

Jo Jo 

for elements Y out of the Banach space 

JJt, p '■= {Y : fl x [0,T] —>• Ei, predictable ||y|| T 2j3 < ooj , T > 0, 
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equipped with the norm defined by 


2 p 

T,2p ■— 


sup E ||Y(t)|| 

0 <t<T L 


2 P 
1 


To show the contraction property of K. on T~Lt,p for small enough T > 0, we first 
decompose 

£(Y)(t) =: S t X o + JC B (Y)(t) + IC c {Y)(t), Y £ H T , P , 


where ICb is the convolution of S with B and /Cc is the stochastic convolution with 
C, respectively. The first term is easiest to handle. From the strong continuity and 
boundedness of S t we get 


(3.13) 


|<S(.)Xo 


I 2 p 
\T,2p 


= sup E 

0<t<T 


||S t X 0 ||f <Ki,iE \\X 0 \\* 


12 p 


where the constant Ki i depends on the bound for the norm of the semigroup. 
The term K, B (Y ) is more difficult to handle. Let Y £ then by Bochner’s 

inequality, Lemma 3.5, Jensen’s inequality and the growth estimate (3.11) 


\\K-b{Y)\\ 2 t \ < sup E 


0 <t<T 


ft \ 2p 

||5 t _ s B(Y(a))|| 1 ds 


(3.14) 


< K sup E 
0 <t<T 




ds 


2p- 


< KT a(l-2p) gup E 
0 <t<T 


(1 + llYWIl!) 


2 P 


/ o 


ds 


< KT 2pa ( 1 + aT~ a sup E 

0 <t<T 


r^)iii 


2 P 


ds 


{t-s) 


1—a 


with constant K changing from line to line, but depending only on p, a and M. Note 
that to apply Jensen’s inequality we have used that (f — s)“ -1 ds is a finite measure 
on (0,t) with mass *“/«, and that the inequality ||a + b\\ 2p < 2 V ^||a|| 2p + ||&|| 2p ^ has 
entered in the last step. Applying the Fubini-Tonelli theorem yields 


sup E 

o <t<T 



< -T a sup E 

OL 0 <t<T 


ll^WIlf 


Inserting into (3.14) we get that 

\\ICb(Y)\\ 2 t p 2p < KT 2pa (l + ||Y||^ 2p ) < oo. 

Let now Y, Z £ TLt, p , then with the same arguments as in (3.14), but with the 
Lipschitz estimate (3.10) instead of (3.11) we obtain 


\\JC B (Y)-ICB(Z)r T p <K! P a sup E 

0<£<T 


\\B(Y(s)) - B(Z(s) 


< A"T (2p - 1)a sup E 

0 <t<T LJo 


(t- s y- a 

f\\Y{s)-Z{s)\\ 

Jo 


- ds 
ds 


2 p 


(t-s) 1 - 


Applying again the Fubini-Tonelli theorem yields 

\\K, b {Y) - IC B (Z)f T p <2p < K"T 2pa \\Y - Z\\ 2p 2p . 
Here, the constants K' and K" depend only on p, a and L. 
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To show similar properties for the stochastic convolution ICc is exactly the same 
as in the proof of the classical result [DPZ14, Theorem 7.2, see pp.189]. Everything 
together yields constants Kr and Kc independent of Xq, such that 

II W) - 1C(Z)\\ Ti2p < II JCb(Y) - Xb(Z )|| Ti2p + II ICc(Y) - X C {Z) || T , 2p 
< {K b + K c )T a \\Y - Z\\ T2p < 1, 

provided T < {Kb + Kc)~ a ■ Hence, K, is a contraction on TLt,p and possesses a 
unique fixed point, which is a mild solution of (3.1) up to time T > 0. Concate¬ 
nating solutions, we obtain a global solution. Finally, to show the the uniqueness 
claim, we consider two arbitrary solutions X\ and X 2 and the stopping times 

tr '■= inf {t > 0 | ||/C s (-X'i)|| 1 > or ^{X^ > R, *e{l,2}}. 

Using the standard procedure as in the proof of [DPZ14, Theorem 7.2], but using 
the estimates for ICb from above and Lemma 3.6 we obtain that the solutions X\ 
and X 2 must coincide up to the stopping time tr. Passing to the limit R —> 00, 
global uniqueness follows. The remaining part, namely showing (3.12) and the 
continuity claim, is subject of Lemma 3.11 and 3.14 below. □ 


Lemma 3.11. Let Assumption 3.1 and (3.11) hold true and let p > 1. Let X be a 
mild solution on [0,T] of (3.1) with initial value Xq £ L 2p {Ll]E\) such that 


E 



\\X{t)\\l P d t 


< 00. 


Then, 


E 


sup ||X(i)|| 

0<t<T 


<K P 


,t( 


1 + E 


2 p 


ll^olli 


If, in addition, St is a contraction semigroup, then the statement is true even for 
P = 1- 


Remark 3.12. We emphasize that the Lipschitz property (3.10) is not needed to 
show this Lemma. 


Proof. We use the notation from the previous proof and write 
X t = S t X 0 + K B {X) t + K c {X) t . 

First, note that the integrability assumption on X and the linear growth prop¬ 
erty (3.11) yield 

(3.15) 


E^ \\B{X{t))\\^ + ||C(^(s)Hhs (! 7 ; £: i) • 

For the case p = 1 we may assume that St is a contraction semigroup. Hence, we 
can apply [DPZ14, Theorem 6.10] which gives 


(3.16) 


E 


r n 


' r T 

sup /CcP0(i)Hi 
,0<£<T 

< KxE 

/ 1 C(^( s )) HS(U.Ei) 

J 0 


For the case p > 1 we use that St is a Co-semigroup and apply [DZ92, Theorem 
1.1] which yields 


r - 


' r T 

sup \\K,c(X){t)\\1 P 

.0 <t<T 

< KpTP-'E 

/ ||C(W(a))||^ s(c/iEi) ds 

J 0 


(3.17) E 

In both cases the growth bound (3.11) yields 


< 00. 


(3.18) E 


IICWWIIhVe,) ds 


< 2 2p M 2p \T 2p + E 


||X( S )||f ds 
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For the drift part we again use the linear growth bound (3.11), and proceeding 
similar to (3.14) we obtain 


sup ||X(s)|| 1 

0 <s<t 

< -S'a.p.T 


2 p 


1 + 11*0 Ilf + 


12 p 


d s 


sup ||-X'(r)|| 1 


2 P 


/O 0<r<s 


- sup I|/Cc(«)|| 1 

0<s<t 


for t < T. Taking expectations, using the Fubini-Tonelli theorem and inserting the 
estimates concerning Kc yields 


E 


sup \\X(a)\\ 

0 <s<t 


<K ( 


t,P,T ( 


1 + E 


ll^olli 


2 V 


E 


sup ||X(r)|| 1 

0 <r<s 


2 p 


ds 


(* —s) 1 


<K a ^ T 1+E ||Xo||i 


II 2p 


E 


E 


sup ||X(s)|| 1 

0 <s<t 


2 V 


ds 


sup ||X(r) 

0<r<s 


Ilf 


ds 

(*-s) 1-c 


Finally, Gronwall’s lemma 3.6 finishes the proof. 


□ 


Ei -continuity of the stochastic convolution K.c(X) t follows from standard results 
and the estimate (3.18). However, for the E\ -continuity of ICb(X) we provide a 
detailed proof, since we can in general not assume that B is Ei-valued. To this 
end, we modify slightly the result [Lun95, Proposition 4.2.1] and its proof. 

Lemma 3.13. Let ip : [0,T] —>• E a be integrable and such that 

sup ||^(s)|| a =: $ < oo. 

0 <t<T 

Then, ( S *ip)t '■= fg St- S ip{s)ds is in C([0, T], E\). 

Proof. Note that for 0 < t < T and arbitrary 0 < e < t, 

-&S t <p = -&S t - e S e <p = ASt-eS e <p = ASt<p, G E, 

where we use that S € ip G T>(A) and f^S t = AS t on V(A); see e.g. [ENOO, Lem. 
II. 1.3]. We now observe that for 0 < s < t < T 

ft fS fS ft — r ft 

/ St-ripir) dr — / Ss-r'ipi'r ) dr = / AS u ip{u) dit dr + / St- r 4’{ r ) dr. 

J 0 J 0 J 0 J s—r Js 

Hence, with Bochner’s inequality and Lemma 3.5 


IKS' *ip) t - (s *^) s ||i 

<K 2>a [ f u a ~ 2 \\^{u)\\ a du dr + K ± . a ( u a ~ x \\^(u)\\ a du 
J 0 J s—r J s 

< f)K a (|t — s|“ + |t“ — s“|) —> 0, as |t — s| —>■ 0. 
For s = 0 we get directly with Bochner inequality and Lemma 3.5 for t\,0, 


S t - r ip{r) dr 


- r i 

< K ha ip (t- r)“ _1 dr = -ipK ha t a -»■ 0. □ 


Lemma 3.14. Under the assumptions of Lemma 3.11 the mild solution X of (3.1) 
is almost surely E\-continuous. 
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Proof. As above, we decompose 

X t = S t X 0 + JC B {X) t + Kc(X) t . 

Continuity of the first summand in the decomposition is immediate, since St is 
strongly continuous. From Lemma 3.11 we get that ip(t) := B{X) t satisfies the 
conditions of Lemma 3.13 and hence continuity of K,B{X) t follows. In the case 
p > 1, estimate (3.18) together with [DZ92, Theorem 1.1] yields continuity of 
Kc{X). For p = 1 we may assume that St is a contraction semigroup on E\ and 
apply [DPZ14, Theorem 6.10] instead. Note that we are always using the continuous 
modifications of the stochastic integrals/convolutions. □ 

Together, these Lemmas complete the proof of Theorem 3.9. 


3.4. Existence of local mild solutions. To obtain only local solutions up to a 
stopping time r, we can relax the assumptions on B and C made in the previous 
subsection. 


Assumption 3.15. There exists a £ (0,1] such that B : E\ —» E a and C : Ei —>• 
HS(t/, Ei) are Lipschitz continuous on bounded sets, i. e. for all N £ N there exists 
Ln such that 

(3.19) || B(Y) - B(Z) || Q + ||C(y) - C(Z)\\ m{U:Ei] < L N ||Y - Z\\ x 
holds for all Y , Z £ E\ with HYH^, \\Z\\ A < (N + 1). 

Remark 3.16. Assumption 3.15 yields for all N £ N constants Mn such that for all 
Y £ Ei with ||F|| 1 <{N + 1) 

(3.20) ||B(Y)|| Q + UCCnilHs^) < M n (1 + \\Y\\ t ). 

Theorem 3.17 (Local Mild Solution of (3.1)). Let Assumption 3.1 and 3.15 hold 
true and let p > 1. Then, for every initial data Xq £ L 2p (PL',Ei) there exists a 
unique mild Ei-solution X of (3.1) on a maximal stochastic interval [0, t[[. More¬ 
over, X is Ei-continuous on [0,t[, t > 0 and lim. t y T ||A(t)|| 1 = oo on {t < oo} 
almost surely. If, in addition, St is a contraction semigroup, then the statement is 
true even for p = 1. 


We use the following localization method, similar to the truncation in [KZS12]. 
For each N £ N fix a monotone decreasing function hw £ with 


(3.21) 


h N (x) 


1, x < N, 

0, x > N + 1, 


and for a constant c > 0, 

sup IIJ^ivL < C. 
N GN 

Define the truncated coefficients 


B n (u) := /iAr(||M|| 1 )S(u), C N (u) := h A r(||u|| 1 )C'(u), 
and consider the localized stochastic evolution equation 
(3.22) dl(")(t) = fil (jv) (f) + B N (X^ N \t))] d t + C N (xW(t)) d W t . 


Lemma 3.18. Let B and C be such that the local Lipschitz assumption 3.15 holds. 
Then, Bn and Cn satisfy the global Lipschitz assumption 3.7 for all N £ N. 
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Proof. First, it is obvious to see that 


(3.23) 


ll-®.v(H)|| a + ||C , Ar(h")|| HS ( [/jSl ) < 


M JV+1 (2 + iV)||y|| 1) 

0, 


<M n (1 + \\Y\\ 1 ). 


||F||i>iV+l 


For the global Lipschitz continuity let Y, Z £ E\ and assume, without loss of 

generality, that ||V|| x > \\Z || r Then, write 

(3.24) 

B n (Y) - B n {Z) = ftivdl^lli) (B(Y) - B(Z)) + B(Z) (Mimii) - Mll^lli)) • 

If ||T^|| x > N + 1, then the first term vanishes. Else, it holds that \\Z 1^ < ||T^|| x < 
N + 1 and thus, in both cases, 


(3.25) IIMnirX^On - B{Z))\\ a < Ln \\Y - Z\\,. 

The second term in (3.24) vanishes if \\Z\\ X > N + 1. Otherwise, 

(3.26) ||B(Z)|| q IMIl^lli) - Mll^lli)l < ^M N {2 + N) \\Y - Z ||,, 

where we applied chain rule and mean value theorem for Frechet derivatives. Of 
course, replacing B by C and ||.|| Q by |MIhs(c/Ei) changes nothing in the com¬ 
putation so that we get a global Lipschitz constant L, depending on c, Mn and 
L n . □ 


For the proof of theorem 3.17 we may assume Xq £ L 2p (fl : Ei), some p > 1, to 
be given and Assumptions 3.1 and 3.15 to be true. For TV £ N we then denote by 
X( N ) the unique mild solution to the localized equation (3.22), which exists due to 
Theorem 3.9. To relax the truncation, we introduce the stopping times 


(3.27) 
and set 


t n := inf {t > 0 X {N \t) ^ > IV } 


(3.28) t := lim tjv- 

N-> oo 

We start with the following preparatory Lemma. 


Lemma 3.19. The stopping times defined in (3.27) and (3.28) have the following 
properties: 

(1) For all k £ N the equality X^ N \t) = X^ N+k ^ (t) holds a.s. for t £ [0, tjv]|. 

(2) The stopping time r is strictly positive. 


Proof. By definition of h n it holds that Bn+i{X^ n \s)) = Bn(X( n \s)) and 
Cn+i(X( n \s)) = Cn{X( n \s)) for s £ [0, tjv]- Hence, using the localization 
property of the stochastic convolution (cf. [BMS05, Appendix A] and [vNVW12, 
Lemma 5.1]), 


X^ N \t) = S t X 0+ f St- s B N (X^ N \s)) d s+ f S t - s C N (X( N \s)) d W{s) 

Jo Jo 

= StXo + f St-sBN+i (AW(s)) d s+ f S t - s C N+1 (X( N \s)) d W(s) 
Jo Jo 


on [0,tjv]. Now, uniqueness of the truncated solutions yields X 1 ,v ) (t) = X( N+1 \t) 
almost surely on [0, tjv]|. For general k £ N the argument can be iterated. 

To show that r is strictly positive, note that it follows from path-wise continuity 
of X( N ) that 


IP [r N > 0] 


lim P 

k—¥ OO 


1' 

tn > r 

= lim P 

sup 

X (JV) (f) 

< N 

k 

k—yoo 

o<t<Vfc 

i 


PIIIAoll, <TV]. 















16 


MARTIN KELLER-RESSEL AND MARVIN S. MULLER 


From the first part of the proof, (tjv) is increasing. Hence, 

P [t > 0] = lim P [t/v > 0] = lim P [||^T 0 IIr < N] = P [|| JSfo || r < oo] = 1, 

N—too N—too 

showing positivity of r. □ 


Proof of Theorem 3.17. For x £ R+ and t £ JO, r[ we set 

(3.29) X(t,x):= lim X^ N \t,x),. 

iV-HX) 

The limit exists, since for almost every ui £ the sequence (X^ N \ui;t,x ))ngn is 
eventually constant for each t £ J0,r[ by Lemma 3.19. It follows immediately that 
1 1—^ X(t,x) is a.s. continuous on each JO, tjv] and hence also continuous on [0, rj. 
Moreover, we may now rewrite tjv as 

t n ■■= inf {t > 0 | ||X(t)|| 1 > N} . 

Continuity of X then implies that the sequence tjv is in fact strictly increasing and 
hence that r is predictable. Moreover, by definition of r we have 

lim ||X(t)|| 1 = lim ||A(tj V )|| 1 = cx) 

t-¥T iV-tCO 

on {r < oo}. 

We focus on the claim that X solves (3.1). By Lemma 3.19 it holds that 
X( N \t) = X(t) on t £ |0,7jv]. Moreover, by construction of B /v and Cjv we 
get B(X(t)) = Bm(X(t)) and C(X(t)) = CN(X(t)) on JO,tjv]- Thus, 

S t X 0 + f St- s B(X(s )) d a + f S t - a C(X)(s)) dW(s) 

Jo Jo 

(3.30) =S t X 0 + f S t - s B N (X( N \s)) ds+ [ S t . a C N (X^(s)) dW(s) 

Jo Jo 

= x (N) (t) 

= x(t) 

holds, and X is a mild solution of (3.1) on [0 ,t/v]. Since N was arbitrary X is a 
mild solution on JO, rj as claimed. To show uniqueness, let X be another local mild 
£i-solution of (3.1) on some stochastic interval JO, rj. For TV £ N we introduce the 
stopping time 


T/v := inf < t > 0 


Xft) > IV or 


X{t) 


> N > At 


Clearly, it holds that 

lim tjv = t A t, a.s. 
N—too 

In addition, for t £ JO, T/v] it holds that 


B(X{t)) = B N {X(t)), and C(X(t)) = C N (X{t)). 


As above, we derive that A is a mild solution of the truncated equation on JO, tjv]. 
The path-wise uniqueness claim of Theorem 3.9 implies X(t) = X(t) for all t, £ 
JO, tjv] and by arbitrariness of N also for t £ J0,t A tJ. Assume now that t < t 
on a set of positive probability. Then lim t/ * T ||X(i)|| 1 = oo on (t < oo} leads to a 
contradiction to the continuity of X. Hence X is unique and [0, tJ is maximal. □ 
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3.5. More global solutions and strong solutions. The results from above can 
be easily extended in two directions: First, we show the existence of global solutions 
under more general conditions, second we show that all mild solutions obtained in 
this section are in fact strong solutions. 

Corollary 3.20 (More global solutions). Let the assumptions of Theorem 3.17 
hold true, but with the local growth condition (3.20) replaced by the global growth 
condition (3.11). Then the solution X is global, i.e. r = oo a.s. 


Proof. By Theorem 3.17 we know that a local solution X to (3.1) exists on a 
maximal stochastic interval [0, r[ and that lim t _j. T ||X(t)|| 1 = oo on {t < oo}. 
Moreover, we know that the stopping time t is the limit of a sequence of stopping 
times tjv < r, that X coincides with the (global) solution X N of the truncated 
equation (3.22) on the stochastic interval [0,tjv] and that ||.X’ jv (tjv) || > N for all 
N £ N. Finally, observe that the coefficients Bn,Cn of the truncated equation 
satisfy the same growth bound as the coefficients of the original equation, i.e. 


(3.31) ||Bjv(y)|| Q + l|Civ(^)|| HS([/iEl) < M (1 + liyy , for all IV eN 

with M independent of N. 

If {r < oo} has measure zero, then r = oo a.s. and the proof is finished. Therefore 
assume, aiming for a contradiction, that P [r < oo] = 2e > 0. By monotone con¬ 
vergence it follows that there exists T > 0 such that P [r < T\ > e from which it 
follows that also P [tn < T] > e for all N gN. Hence 


(3.32) 


E 


sup \\X N (t)\\ 2 1 P 

0<t<T 


> N 2p V [t n <T}> N 2p e 


holds. On the other hand, applying the growth bound (3.31) and Lemma 3.11 to 
each Xjy it follows that 


E 


sup \\X N (t)\\ 2 * 
0 <t<T 


< K p ^t ^1 + E 



with the right hand side independent of N. Combining with (3.32) and choosing 
N large enough, the desired contradiction is obtained. □ 


Corollary 3.21 (Strong solution). Under the assumptions of Theorem 3.9 and 
Theorem 3.17 respectively, the processes X^ N \ N £ N and X are also the unique 
strong solution of respectively (3.22) and (3.1). 

Proof. Let (X,t) be the unique mild solution from Theorem 3.17, and X^ N \ tn 
and respectively Bn and Cfy be the solution, stopping times and parameters of the 
truncated equation (3.22), corresponding to (3.1). By the natural embedding we can 
identify the Ify-paths of with paths in E. Further, Kuratowski’s theorem (cf. 
[Kur66]) implies for the corresponding Borel cr-algebras that © (E±) = © (E) D Ei 
and hence that we can extend Bn trivially to a Borel function on E without affecting 
the regularity properties that Bn has on E±. Writing down the Hilbert-Schmidt 
norm one immediately observes that also 

HS(C/,.Ei) HS (U,E). 

Both are seperable Hilbert spaces so that we can argue in the same way to extend 
Cn to a Borel function from E into HS(f7, E). The stochastic evolution equations 
now fit in the framework of [PR07, Appendix F], where sufficient conditions for 
obtaining weak and strong solutions from mild ones are given. Proving the corollary 
now simply amounts to showing that these conditions are satisfied. 
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First, recall (3.15) which yields for all T > 0 that 


[ B n (X^ n \s) + C N (XW( S ) 
Jo 


HS (U;E) 


ds < CXD, 


-almost surely and from equivalence of the norms of E\ and T>(A) we get a. s. 


f X( N \s) + AXW(s) 
Jo 


ds < CXD. 


For the step from mild to weak solutions we also need to verify for all g £ V(A*) 


(3.33) 


E 


S t C N (X^ N \s)),A*g 


HS(C7,R) 


ds df < CXD. 


Recall that A generates a strongly continuous semigroup on E. Using Cauchy- 
Schwarz inequality we then get for any complete orthonormal system (e*,) of U and 
0 < s < t 


(s t - B C N (xW(s)),A*g) E 


2 

HS(t/,R) 


J2\{St-sC N (X( N \s))e k ,A*g) E 


< 


< 


S t - a C N (xW(8)) 


HS (U;E) 


\\A*g\\ 


C N (xW(s)) 


HS (U,E) 


\\A*g || 


E ' 


2 

E 


From (3.15) (or boundedness of Cjv) it follows that (3.33) indeed holds true. As 
we have seen in the proof of the continuity part of Theorem 3.9, 


[ S t - s B N (X( N \s))ds, and t^f S t - s C N {X^ N \s)) dW s 
Jo Jo 

are continuous, and therefore predictable. Hence, all integrability and measurability 
assumptions which are needed to apply [PR07, Prop. F.0.4 and F.0.5] are satisfied 
and we conclude that the mild solution X^ is also a weak and strong solution 
of the truncated equation (3.22). Applying the the same localization argument as 
above, but now on E, the result translates to X. 

For the uniqueness claim, suppose first that is another global strong solu¬ 

tion of the truncated equation (3.22). Applying the results in [PR07, Appendix F] 
we get that Y < - N ' 1 is also an .E-mild solution. Since Bn and Cn are bounded, we 
get that 



S t - s B N (Y^)(s) 


Ei 


S t - s C N (Y^ N \s)) 


2 

HS(U;Ei) 


ds < CXD 


almost surely, so that yW is even an £d-mild solution. Hence, the uniqueness 
part of Theorem 3.9 yields yW = X^ N ^ almost surely. Finally, for another local 
strong solution (Y. <,'), truncating with respect to the El-norm yields Y = X on 
PUAt[, almost surely. Since [0,r[[ is maximal for the Ei-mild solution, we obtain 
that <d < t almost surely. □ 


4. The fixed BOUNDARY problem as a stochastic evolution equation 

The goal of this section is to reformulate the free boundary problem (2.1) as the 
abstract evolution equation (3.1). We start start by identifying the appropriate 
function spaces and introduce 

£ 2 := T 2 (K+) © L 2 (R + ) 0 K, 


:= H k (R+) © H k ( R+) 0 R, 
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where, as usual, L 2 denotes the Lebesgue space, H k the fc-th order Sobolev space. 
Recall that © denotes the direct sum of Hilbert spaces, i.e. the scalar product on 
£ 2 is defined through the scalar product on L 2 (R+) by 

{(u 1 ,u 2 ,x),(v 1 ,v 2 ,y )) s2 = (Mi,^i)i2 (R+) + (u 2 ,v 2 ) L 2 (M+) +xy, 

and similarly for the space S) k . 

Recall the definitions of the operators A, B , C that were given in equations (2.11), 
(2.12), (2.13) in terms of p±, a±, p and </. To define the domain T>(A) of A we set 

D := H 2 (WL+) n H](R + ), and V(A) =DxflxICfi 2 . 

Finally, the space 'D(A) shall be equipped with the graph norm 

IML := IML* + ll-'MU*. u e v{A). 

Note that on T>(A), the graph norm is equivalent to the Jo 2 -norm, as can be seen 
from integration by parts and the Cauchy-Schwarz inequality. Moreover, T>(A) is 
a closed subset of Sj 2 . 


To provide the connection with the results of section 3, we set 
E 0 = 2 2 , 

As the norms of T>(A) and Sj 2 are equivalent we may use either one to topologize 
E i- The following result holds true for interpolation between T>(A) and £ 2 : 

Lemma 4.1. For a £ (0,1/4) it holds that, 

E a = £) 2a , 

with equivalent norms. 

This follows from the fact that Z?((c — A)“) = U 2 “(R + ) for the Dirichlet Lapla- 
cian A, iff a < 1 /a. To our knowledge, this was shown first in [Gri66], but see 
also [LM72, Ch.l , Thm. 11.6]. By the structure of V(A) and £ 2 , this directly lifts 
to E a . This result can be understood in the sense that the boundary conditions, 
which distinguish T>(A) from f) 2 ‘are lost’ during interpolation exactly at a = 1/4. 


Lemma 4.2. The operator A, defined in (2.11), satisfies Assumption 3.1. 


Proof. The Dirichlet Laplacian on L 2 (R + ) is a self-adjoint operator. This property 
is inherited by A , which is hence a self-adjoint operator . Moreover, 


- (Au,u) £2 = ||^ui| 


L 2 ( R+) 


+ ll ^ 2 " 2 


L 2 (R+) 


c\\u\\l 2 >c\\u\\% 


for all u £ T>(A) and hence —A is positive in the sense of [Lun09]. By [Lun09, 
Lem. 4.31] it follows that (—c, oo) is contained in the resolvent set of A and the 
resolvent satisfies 

||R(A,A)|| £2 <-^, a > —c. 

Choosing M > max(l, 1/c), the estimate (3.2) follows and Assumption 3.1 is satis¬ 
fied. □ 


Lemma 4.3. Suppose that p±, a±, p and £ satisfy Assumptions 2. 2, 2. 3, 2.5 and 
2.6. Then the operators B(u ) and C(u), defined in (2.12) and (2.13) satisfy the 
local Lipschitz Assumption 3.15. 

Proof. We decompose B{u) into 

(4.1) B(u) := S M (u) + B p (u) + cu, uGV(A), 
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where 

/ M+(-) u i> ^ u i) \ / \ 

B p {u) = I ^U 2 ) I , B p {u) = p{X{u))Xu and Vm := I -^ 1 x 2 I 

Recall from (2.10) that I (it) is the vector of boundary values given by I{u) = 

0+), ^rt 2 (f, 0+)). The trace operator is known to be continuous on H 2 { R + ), 
see [LM72], and hence, by equivalence of norms we have that I £ L(ij 2 ,R 2 ) with 
operator norm Kx, say. Denote by .B/v+i the (closed) ball of radius N + 1 in Sf 2 . 

The image of -Bat+i under X is closed and bounded, hence a compact subset of R 2 . 

By Assumption 2.5, the function p : R 2 i-a R is locally Lipschitz continuous and 
hence Lipschitz on any compact subset of R 2 . Thus, we find a constant L p ,a such 
that 

\p(l(u)) - p(T(w))\ < L PiN \X(u) -I{w) I = L PiN K x \\u-w\\^ , Va,iu £ Bjv+l 

By a similar argument and using only continuity of p instead of the Lipschitz 
property, we find M p> jv such that 

sup \p(X(u))\ < M PtN K x . 
il£Bn+i 

Finally, for any u,w £ -Bjv+i and setting Lm = Kx{M p ,n + {N + 1)L Pi jv) we 
obtain 

|| p(l(u))Vu - p(I(w))Vw|| Ql < \p(T(u))\ || Vu - Vw|| ftl + ||Vu|| ftl \p(X(u)) - p(X(w))\ < 

< M P)N K x ||ix - w\\^ 2 + {N + 1 )L p , n Ki ||ix - w\\^ 2 = 

= L N \\u-w \\^ 2 ■ 

Now let a < 1/4. By the continuous embedding fj 1 f) 2 “ and the equivalence of 
norms on 55 2 and V(A) we may find K a such that also 

|| p(I{u))S7u - p(l(w))\?w ||^ 2c< < K a L N ||ix - . 

Recalling that fj 2a = E a and X>(A) = E\, this yields that B p is Lipschitz as a 
mapping from bounded subsets of E\ to E a . 

It remains to show the same properties for B p and for C, which we delay to 
Appendix A and B, respectively. □ 

In particular, with the results from section 3 we get a unique solution X = 
( 111 , 1 x 2 , 2 ;*) of (2.9) on the maximal interval [0,r[. Applying the following lemma 
to Corollary 3.20 yields global existence, under global growth assumptions. 

Lemma 4.4. Suppose that in addition to the assumptions of Lemma f.3 also As¬ 
sumption 2.13 holds and p is bounded. Then the operators B and C satisfy the global 
growth bound (3.11). 

Proof. Decompose B into B p and B p as in the proof of Lemma 4.3. Using Assump¬ 
tion 2.13 we get the point-wise estimate 

\B p (u)i\ < (a + b(ui, -§^Ui))(\ui\ + \&Ui\)) , i £ {1,2} 
where a £ L 2 (R + ) and b £ L°°(R). Taking £ 2 -norms yields 
(4-2) \\B p (u )\\ s2 < 2 (Halloa + H^Hoo IMIflO • 

For the first weak derivative we extract from the proof of Theorem A.6 (cf. Eq. (A.l).) 
-£h+(x,u i(x), ^ui(x)) =-§ i p + {x,u 1 (x), -§^ui(x)) 

+ JEui(x)) (a;) 

+ -§zP+{x,ui{x), -j^ui(x))-j^ 2 U\{x), x £ R+. 


(4.3) 


A STEFAN-TYPE STOCHASTIC MOVING BOUNDARY PROBLEM 


21 


Whereas the first summand admits a bound similar to (4.2), part (b) of Assump¬ 
tion 2.13 yields for some b 6 L°°(l R 2 ;R), 


^// + (x,tii(x), J^i(x))£m(x) < b(ui(x), £u 2 (x)) \£m 


and 


d 


a^lJ-+(x,u i(x), ^ui(x))^ui(x) < 6(ui(x), ^u 2 (x)) 

Of course, the same holds for p_ and u 2 so that we can summarize 

b 


d 


(4.4) 


d 

dtc 


BM\ 


02 <2 


li 2 


Collecting all the estimates and using S) 2 c —> f) 1 we get a constant M^, depending 
on a, b and b only, such that 

(4-5) WB^W^KM^l+WuW^). 

Using the assumption that p is bounded by a constant M p we easily estimate 
II^WII^i = ||p(I(u))Vu || J - 1 < M p [| U || fl2 . 

Combining with (4.5) we obtain 

||B(«)|| fll <(Af M + M p + c)(H-W fla ). 

Using the continuous embedding A 1 Sj 2a and the equivalence of norms on H 2 
and V(A) as in the proof of Lemma 4.3 yields the global growth bound (3.11) for 
B. 

To show the analogous growth bound for C, observe that the following equalities 
hold for all x £ R + and u £ 13(A) C f) 2 : 

£a+(x,u i(x)) = £a]_(x) + £al(x)u!(x) + al(x)£u!(x) 

£z<j + (x,u i(x)) = £zo\{x) + -£?<tI(x)u i(x) 

+ 2 m a K x )m u ^x) + o-l(x)£zu i(s). 

Hence, there exists a constant K > 0 such that 

ll f7 +0) U 1 (•)) llff2(R + ) < || °"+ l|pf 2 (R + ) 4” ^ || <7 + || C 2 (R + ) ll-ff 2 (R+) ‘ 

We apply the same argument to er_(—., it 2 (.)) to obtain 


(4.6) 


/ cr + (.,Ui(.)) \ 

I a-(-.,u 2 (.)) I 


— Ka (1 + IMI^) , 

A 


for a constant K a , depending on and <r_ only. By Assumption 2.3 tJ±(.,rti(.)) 
satisfies Dirichlet boundary conditions at 0 and we may apply Lemma B.4 to obtain 


II^( M )IIhS(C/;X>(A)) — ^ 


sup 

xgR 


E 

i— 1 


k (i) ow) 


L 2 


/ a + (.,u i(.))\ 

( a-(-.,u 2 (.)) j 


A 


which together with (4.6) yields the desired linear growth bound. 


□ 


Finally, we prove a refined result on the blow-up behavior of the solution X(t) 
at the stopping time r. We show that under Assumption 2.13 a finite-time blow-up 
of X(t) can only happen if the boundary values l(X(t )) themselves blow up. For 
IVgNwe introduce the stopping times 

To ,at := inf {t > 0 1 1 < t, \l(X(t))\ > N} . 

and set 

r 0 = limsupr 0) Ar. 
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Recall at this point the convention that inf 0 = oo. From Theorem 3.17 we know 
that 

(4.7) lim ||X(T)|| . = oo, a. s. on (r < oo}. 

t/'T 

Since a blow up of \I(X(t))\ implies a blow-up of the norm ||X(t)|| w2 , and hence 
also of ||X(t )||_4 we obtain that r < r 0 a.s. and hence only two events are possible: 
Either 

> r = r 0 , i.e. a blow-up of X(t) coincides with the blow up of the boundary 
values I(X(t)), or 

> t < oo, but t 0 = Too, i.e. a blow-up of X(t) occurs without simultaneous 
blow-up of its boundary values. 

The following theorem shows that Assumption 2.13 rules out the second case: 


Theorem 4.5. If in addition to the assumptions of Lemma f.3 also Assump¬ 
tion 2.13 holds, then 

P[ To =t]=1. 


Proof. Because of r 0 > t and maximality of r, (4.7), it suffices to show 
(4.8) lim ||X(t)||. < oo, on {r < oo}. 

t/'To.N^T 


Indeed, this yields t 0j n < x on {r < oo} and thus r 0 < r almost surely. For N £ N, 
let Ln : R —> R be the truncation function defined in (3.21) and define 

PN{x,y) := p(x,y)h N (\(x,y)\), (a :,y) G R 2 , TV G N. 


Then, due to Theorem 2.14, equation (2.9) with p replaced by p^ admits a unique 
global solution denoted by Xjj. Since p = p^ on the ball of radius TV, we have 
X{t) = TVat(T) on JO, r 0 ,iv AtJ. Since Xn is D(A)-continuous, we have 


lim ||X(i)||^= lim ||X A r(T)||_ A < oo, 

£/*T 0 ,Ar At t/^To^ At 


on {t < oo} 


for all TV G N and the proof is complete. 


□ 


5. Transformation from moving to fixed boundary 

As the last step towards a complete proof of the main result Theorem 2.11 we 
make the transformation (2.7) to the fixed-boundary equation (2.8) rigorous. Since 
the equation is stochastic and its solution not differentiable in time, the classic 
chain rule cannot be applied. As an alternative we could use Ito’s formula, which 
requires the transformation map to be C 2 . It turns out that the transformation is 
only C 1 , but linearity in its first argument and the bounded variation of x* (t) are 
in combination sufficient to make a stochastic version of the chain rule work. 


5.1. Stochastic chain rule. For a given cylindrical Wiener process W on U we 
consider the H 1 (R)®R-continuous process {v, x) on JO, t[, for a predictable stopping 
time t, such that, 


(5.1) 


V t =v 0 + p. s ds+ a s dW s 

Jo Jo 

r t 


s(t) =x 0 + x s ds, 
Jo 


where /i : [0, t[— >• L 2 (R) and a : |0,t[— >• HS(f7, L 2 (R)) are predictable processes, 
and x : JO,t[— > R is continuous and adapted. 

For iGlwe define the shift operator 9 X acting on L 2 ( R) as 


(5.2) 9 X : L 2 (R) —» L 2 (R); /(.)->/(a: + .). 
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Observe that the shift operator is a linear isometry and hence continuous. It is 
obvious that the shift operators (0 x )xeR form a group under composition, i.e. 9 X 9 f = 
9 X+ ^; in fact this group is strongly continuous in L 2 (R), in the sense that 

(5.3) lim || 9 x f - /|| i2 = 0, for all / £ L 2 (R), 

x^f-0 

see e. g. [Wer07, Section VII.4]. The same properties hold true for the restriction 
of (^x)xeR to the Sobolev space J7 1 (M). Finally, consider the function 

(5.4) F : L 2 (K + ) © R —>• L 2 (R), {v,x)^9 x v 

which formally transforms the solution (u,x*) of the fixed boundary problem (2.8) 
into the solution (F(v,x*),x*) of the moving boundary problem (2.1). 

We start with a Lemma on some uniform continuity estimates for the shift op¬ 
erator on L 2 (R). 


Lemma 5.1. Let T > 0, s i-A f s and s i-A g s be continuous functions from [0,T] 
to L 2 ( R) and HS(Z7, L 2 (R)) respectively, and let 9 X be the shift operator on L 2 (R). 
Then the following holds: 

(a) For every e > 0 there exists S > 0 such that 

\\9hfs — ft\\L 2 (K.) — e 

for all \h\ < S and s,t £ [0, T] with |s —1| < S. 

(b) For every e > 0 there exists 5 > 0 such that 

II 9h,g s — 9t IIhs(u,l 2 (r)) — e 
for all \h\ < S and s,t £ [0, T] with |s — t| < 6. 


Remark 5.2. Applying the lemma to the constant function g t = g, 9 can be consid¬ 
ered as a strongly continuous group on HS(C7; L 2 (R)). 

Proof. For claim (a), note that due to the continuity of s K>• f s there exists N := 
N(e) £ N and t \,..., tjv G [0, T\ such that 

■ /f in vr ~ - I Vs e 


i.e., we can find N balls in L 2 (R) around the points f ti of radius which cover 
the whole range of f s . Moreover, due to the strong continuity of the group (0 x ) x gR 
there exists S > 0 such that 


\\0 h fu-ft J La( R)<! V\h\<5,ie{l,...,N}. 

The estimate 

\\9hfs — /t llz, 2 (M) — — /*« IIl 2 (R) + ll/s - fti IIl 2 (R) + II ft ~ /*«IIl 2 (R)} 



concludes the proof of the first claim. 

We show (b) in a similar way; to alleviate notation we denote by ||.|| ffjS . the Hilbert- 
Schmidt norm ||-|Ihs(I 7 L 2 (r))- D ue f° th e continuity of s i->- g s there exists N := 
N[e) £ N and t \,..., tjv G [0, T] such that 


min 

ie{i,...,iv} 


Us ® 9u\\hs 



Vs £ [0, T). 


Denote by (ek)keN an arbitrary orthonormal basis of the (separable) Hilbert space 
U. Since g s is a Hilbert-Schmidt-operator for every s £ [0,T], the series ILgtJljys = 
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llffu e fe|lL 2 (R) converges for every , N}. Hence, there is M := M(e) G N 

such that 

oo 2 

X! \\9uek\\ 2 L 2 im < — Vi G {1,..., N}. 

k—M+l 

Moreover, due to the strong continuity of the group (0 x ) x gR there exists S > 0 such 
that 

m e 2 

Y II e h(gue k ) - 9ue k \\ 2 L 2 m < — V\h\ < 6, i G {1,..., N}. 

k =1 

Combining with the previous equation, we obtain 


II ®h9ti dtiWnS ~~ v II QhjgtjCk) #£i e fc|lL 2 (R) < 


k =1 
M 


— 5^11 OhifftiCk) #£i e fc||.L 2 (R) + ^ ll^*< e fc|lL 2 (R) — 

k—M-\-l 


k=1 


S I2 +4 38 = T Vi 6 {1,... ,N}. 


The estimate 


\\0h9s 9t\\ns — /P* 11 , n { dtiWns Wds Su IIlts + lll?t ^JIlts} 


<3| — + -—h — ) =e z 

“ 6 12 12 , 


concludes the proof. 


< 


□ 


For notational simplicity, we denote by v' the weak derivative of v G H 1 ( R). 

Lemma 5.3. The transformation F from (5.4), restricted to H 1 { R) © R, has the 
following properties: 

(1) F is a continuous mapping from H l { R) © R to i? 1 (R); 

(2) F is a continuously differentiable mapping from iJ^R) © R to L 2 (R) with 
Frechet derivative given by 

(5.5) D {VtX) F{h,£)=0 x h + ti6 x v'. 

for i,(£R and v, h € FI 1 (R). 

Proof. For (1), we estimate 

\\0 x u — ^ |^ X ( U — ||#1(r) = ll w — x w ll.£fi(R) — 

— II M — 17 II AT 1 (1R) + \\ V ~ ^£-x^||jji(r) ) 

where by the strong continuity of 9 X on F[ 1 ( R) the latter term vanishes as |£— x\ —> 0. 
To show (2), we first verify that (5.5) gives the Frechet derivate of F, by estimating 

R 2 ■= ||(F(u + h,x + f) - F(v,x) - D {VjX) F(^,h) ||* 2(R) = 

— 11 ®x-\-£ “1“ £,@xV ^ :c ^'IIl 2 (R) — 

— 2 \\@€ v ~ v ~ £ v IIl 2 (r) ^ l|0f^ — ^IIl 2 (r) = 

poo poo 

= 2 / (v(y + 0 - v(y) - £v'(y)) 2 dy + 2 ( h(y + £) - h(y)f d y. 

9—00 


—oo 
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Applying first the fundamental theorem of calculus and in the second step Jensen’s 
inequality and Fubini’s theorem we continue with 


R 2 < 2£ 


I (/ ^ V '( y + ~ V '^ d ^) d?/ + / (/ h '( y + z ^ dz ' S ) dy 

— ^ ll^2f u ~ v IIl 2 (R) + 11^' lli 2 (R)| — 


< 


< 0 (l£| 2 ) + 2£ 2 


showing Frechet differentiability of F. To show continuous differentiability we es¬ 
timate 


||D(v,x)F(£, h) D w ^ y F^, h )||— ||£$x^ d~ 9 x h £,9 y w y yh\\L 2 (R) — 

A |£| || O x —yV — W ||£2( R ) + || Ox—yh — h \\£ 2 ( R ) A 

< |£| \\e x - y v' -w'\\ L2{m + \x - y\ ||/i'|| i2(R) < 

< (l£l + ll^llffi(R)) max (ll Ox-yv' - w'|| i2(R) ,\x-y 
Writing ||.|| for the operator norm from iJ 1 (M) ©M to L 2 (R), this shows that 

||-D( l)> x)-P 1 — D(w,y)F\\op — max (\\9 x — y V ~ w ||£ 2 ( R ) , \x — y |^ . 

The right hand side goes to zero as ||w — w|l.ffi(R) + \x — y\ -A 0 by Lemma 5.1, and 
hence D^ V X ^F depends continuously on ( v,x ). □ 


Theorem 5.4 (Stochastic Chain Rule). Let (v,x) be given by (5.1) and set Ut = 
F(vt,Xt) = 0 Xt Vt . Then u satisfies 

(.(fi Xa IJ>s) + x s (9 Xs v' s )) d s+ f (9 Xs a s )dW s , 

Jo 

on [0,t[, where the first integral is an L 2 -Bochner integral and the second one an 
HS -stochastic integral. 


(5.6) 


u t = u 0 + 


Proof. Without loss of generality assume r = oo almost surely. Else, take an 
announcing sequence of stopping times (tjv) for r. Then, multiply x, yt and at 
with the indicator function 1 jo, TJV j and replace resp. a;* and v by x,/\ TN and u Arjv . 
In this proof we use the shorthand notation x s< t = Xt — x s , v s j = vt — v s etc. For 
V an arbitrary partition of [0,T] we decompose 

F(v t ,x t )-F(v 0 ,x 0 )= ^2 F(v t ,x t ) ~ F(v s ,x s ) = 

[s,t]eP 

= ^2 D (x a ,v s ) F (Vs,t,X s ,t) + R{v s ,X s ,V s>t ,X s>t ). 

[s,t]ev 


where DF is the Frechet derivative of F from (5.5) and R is a remainder term. 
Using equation (5.1) and the explicit forms of F and DF we rewrite 


D (v s ,x s } F s,ti x s ^t) — 9 Xs Vs,t + x s ^t9 Xs v s — 

= [ (Qx a Vr) dr+ f x r (0 Xs v' s ) dr + f ( 9 Xs a r ) dW r 
J s J s J s 

R(v S ) 'Esi Vs,t’> — ^Xt^t ^Xs^t %S,t (^Xg^s) • 
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The formula (5.6) follows, if we can show each of the following (a.s.) convergence 
statements 


(5.7) 

(5.8) 

(5.9) 


and 

(5.10) 


E / ^ (^X) dr 

,t]ev n Js 

T. / ( 0 Xa fji r ) dr 

,1 ^rr~, J S 


s,t]ep n 


/ X r ( 9 Xr v ' r ) dr 
Jo 


(0 Xr Hr) dr 


^ ) {0 Xt Vt 9 Xs V s X s ^t (^i s ^ s )) t 0 
[s,t]e-p„ 


E 

s.tje-p, 


f (0 Xs cr r ) dW r [ ( 9 Xr a r ) dW r 

Js Jo 


for some sequence of partitions {V n )ne n as n —> oo. 

For claim (5.7) chose e' > 0 and set e = e'/14[0,T] where 14[0,T] is the total 
variation of the process x over the interval [0, T], By Lemma 5.1 we can find 5 > 0 
such that 

\\8hV s — II i 2 (M) — e 

for all \h\ < 5, |s — r| < S. By continuity of x we can choose the mesh of V n fine 
enough such that |x rjS | <8 for |s — r| < \V n \ and hence 


[s,t]eP, 


/ * 


Xr (0x a v' s - 0 Xr v'r) dr 


< 


L 2 ( R) 


< 


E 

[s,t\GV n 


.V s - v r 


Il 2 (r) 


dr < V x [0 ,T]e = 


for \V n \ small enough. As e' > 0 was arbitrary this shows (5.7) for any sequence of 
partitions with \P n \ —>• 0. For (5.8) we estimate 


E /' 

[s,t]GV n Js 


{@x a Hr 


0x r Hr) dr 


< 


\ ' 


L 2 (R) [s.de'Pn 


II ^X a —Xr\^T 


f^ r IIl 2 (r) dr. 


The integrand converges to 0 for \V n \ \ 0, due to strong continuity of 9 and 
continuity of x. By dominated convergence, this carries over to the whole integral. 
For (5.9) estimate 


R := 

E 

( @Xt^t Qxs'Vt %S,t (^Cs^s)) 


[s,t]ev. 


< 

E 1 

^X S) t V t — Vt — Xg^Vg II = 




= 

E ( 

/ r°o 

/ (vt(y + x 8i t) -v t (y) -X 

I s 

,t]eP„ 

W — oo 


< 


m 


1/2 
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Applying first the fundamental theorem of calculus and in the second step Jensen’s 
inequality and Fubini’s theorem we continue with 

R = E \ x *a([ [ ( v 't(y + zx s,t) - v'siy)) 2 dy dz\ < 

+1 (zHD \j —oo j 0 / 


[s,t]ev, 

< E 

[s,t]ev n 


yt-<\\ L 2 (WL) & 


1/2 


Again yh p |x s ,f | can be bounded by the total variation of x, while Lemma 5.1 
shows that the L 2 -norm vanishes uniformly as \V n \ —>• 0. Finally, to obtain the 
convergence of the stochastic integrals in (5.10), we define the HS(I7, L 2 (R))-valued 
functions 

^P,W= E ( e xa a r) l(M]( r )> ®{r) = 0 Xr o r 

[s,t]eV n 

and note that (5.10) is equivalent to lim^oo 4>-p n (r) d W r = / Q T $(r) d W r along 
the sequence of partitions V n . Denoting by ||.|| ffS the Hilbert-Schmidt-norm on 
HS(t/,L 2 (M)) we claim that 


(5.11) 



||<Fp„(r) - <f>(r)f HS 


dr —>• 0, 


Indeed, by rewriting 


as [P n \ —> 0. 



II ®v n {r) - $(r)llL? 


dr = 


E 

[s.iJe'Pn 


|| 0x a -x 


(J y 


2 

HS 


dr 


we may proceed as in the case of (5.8), since 9 is strongly continuous on HS([/; L 2 (R)), 
see Remark 5.2, to conclude that the Hilbert-Schmidt-norm on the right hand side 
vanishes uniformly as \V n \ —>• 0. Choosing <5, e > 0, [DPZ14, Prop. 4.31] says that 


/ 


d ) 'Pn( r ) ~~ $(»") d W r 


>6 


L 2 (R) 

< e + ■ 


^ ll^Pn M - $ WllE dr> ^2 


and we see that as \P n \ —» 0 the right hand side can be made arbitrarily small. 
Hence, 4>-p n (r) dW r f () 4>(r) dW r in probability, as || —> 0. In particular any 
sequence of partitions with mesh tending to zero contains a subsequence ('Pn fc )fceN 
such that the convergence of stochastic integrals takes place almost surely, com¬ 
pleting the proof of (5.10). □ 


5.2. Completing the proof of the main result. Making use of the notation 
defined in section 2 and 4 we complete the proof of Theorem 2.11 by combining the 
relevant results. 


Proof of Theorem 2.11. Given initial data xq £ M, vo £ T(a;o) we set 

X 0 := (v 0 (x 0 + (.))|r> 0 ,vo(x 0 - (-))k> 0 >^o) £ T>(A). 

By application of the Lemmas in section 4 to Theorem 3.17 and Corollary 3.21 
there exists a unique maximal strong solution of (2.9) X on £ 2 with initial value 
Xq. Note that X is f) 2 continuous and denote its .^-explosion time by r. Setting 

ui(t) := Xx(t), u 2 {t) := X 2 (t), x*(t) := X 3 (t) 
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we obtain a solution of the fixed boundary problem (2.8) on |0,t[. Next we paste 
together u\ and u 2 by setting 

u t {x) = ui(t, x)1r + (x) + u 2 (t, -x)lu_ ( x ) 

and note that due to the Dirichlet boundary condition at 0 it holds that Ut G 
id 1 (R). Thus we may apply the stochastic chain rule of Theorem 5.4 to vt '■= 
F(ut, —x*(t)) = 0— x *(t)Ut t° obtain a local solution v t . of the stochastic free bound¬ 
ary problem (2.1). 

To show uniqueness we assume that there exists another local solution (v,x*) of 
(2.1) on a stochastic interval [0,c[, which, by definition is a. s. H A (M)-continuous. 
Applying the stochastic chain rule of Theorem 5.4 to Ut = F(v t , x*(i)) = Ox,(t)Vt we 
obtain a local solution u t of the stochastic fixed boundary problem (2.1). Reversing 
the procedure from above ut can be rewritten as a solution X(t) of the abstract 
stochastic evolution equation (2.9). The parts of Theorem 3.17 and Corollary 3.21 
on uniqueness and maximality imply ? < r and that X[t) is equal to the solution 
X(t) constructed previously. We conclude that also (v,x*) = (v,x*) and the proof 
is complete. □ 


Appendix A. The Nemytskii operator on Sobolev spaces 


In this section we prove some regularity results on the Nemytskii operator N, 

N(u)(x ) = p(x,u(x)), 

on the Sobolev spaces H k ( R+). Here, p : R + x —s- I and x £ R+. Note that 

these results are well-known, even for more general spaces, in the case of bounded 
domains, see e.g. [Val88], [AZ90]. However, in the case of unbounded domains 
several additional conditions on /i are necessary to make them work. First, we 
state a result which guarantees, that under certain assumptions on /x, N maps H k 
into H k . For a proof we refer to [Val85, Theorem 1], of which it is a special case. 

Lemma A.l. For each integer k > 1 the space H k { R + ) is a Banach algebra. In 
particular, there exists a constant c such that for all u, v G H k (WL + ) it holds that 
uv G H k { K + ) and 

\\uv\\ Hk <c\\u\\ Hk . 


Next, we adapt [Val85, Theorem 2] to our setting. For notational reasons we 
also introduce the Nemytskii operators 

N x (u)(x) := (x,u(x)), N y .(u)(x) := ( x ’ u ( x ))> 3 = 

for u G H k ( R + ; R d ), x G R+. In order for N to map H k into H k again, we need 
certain growths restrictions, which is not the case on bounded domains. 


Assumption A.2. Assume p G C' m (R>o x R d ,K) and 

(a) For each integer l, 0 < l < m there exists an ai G L 2 (R+) and some 
bi : R d —>• R+ locally bounded, such that 


p(x,y) < ai(x) + bi(y) \y\ , ViGR + , 1 / 


(b) For each multiindex a with or ^ |a| < m, the functions sup^-gg \D a p(x, .)| 
are locally bounded. 


Assumption A.3. Assume that p G C rn (M>o x M d ,R) and D a p(x,.) is locally 
Lipschitz for all multi-indices a, |a| < m with Lipschitz constants uniform in x G 
K>o, i- e. we assume that for all r > 0 there exists L r > 0 such that 

I D a p(x, y ) - D a p(x, z)\<L r \y-z\. 
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holds for all y, z E R d with \y\, \z\ < r and a, |a| < m. 

Remark A.4. If y satisfies Assumption A.2 for some integer m > 1, then y satisfies 
Assumption A.3 for m — 1. 

Remark A.5. Recall the Sobolev embeddings 

H m+1 ( R+) -A BUC m ( R+), 

where BUC m (R + ) denotes the Banach space of functions with bounded and uni¬ 
formly continuous derivatives up to order m. As usual BUC m ( M+) is equipped 
with the C' m -norm. In the following, we will work with the BUC m representative 
of the elements in H m+1 without further comment. 

Theorem A.6. If Assumption A.2 holds for some integer m > 1, then the operator 
N is continuous from |_)) d into H m ( R+). 

Proof We adapt the proof of [Val85, Theorem 2] for the domain R + and the spaces 
H k by incorporating the additional growths assumptions. We proceed by induction 
and consider m = 1, first. Since y E C 1 (M>o x R d ) we get immediately that 
N(u), N x (u) and N Vj (u ), j = 1 are bounded and continuous functions for 
u E H 1 (R+,R d ) fixed. Let now ( u n ) C 7L 1 (R+; R d ) nC 00 ^;^) such that 
u n — > u in H 1 . Then the convergence also takes place in ||.|| and by the chain 
rule we can write 

n 

(A.l) ^(« n ) = N x {u n ) +Y, N yj (u n )Vu™. 

i =i 

By assumption, 

|A r (w n )| < a 0 + (b 0 (u n )) \u n \, and |iV I (u")| < 0l + (6i(«")) \u n \ . 

Since u n —> u uniformly and in L 2 we get for each estimate that ||&i(it ra )|j for 
i E {0,1} is bounded in n € N such that both, N(u n ) and N x (u n ) are bounded by 
T 2 -converging sequences. Hence, we can apply a version of Lebesgue’s dominated 
convergence [Kal02, Theorem 1.21] and obtain the L 2 convergence of 

N(u n ) —> N(u), and N x (u n ) —» N x (u ), as n —>• oo. 

For the remaining summands we get 

|| N yj (u n )Vu] - N yj (u)Vu\\ L2 < sup ||iV % (u n )|| oo \\Vu] - V Uj \\ L2 

n£ N 

+ 11 {N Vj (u n )-N yj (u)) Vu 0 \\ L2 

which goes to 0 as n — > oo. Indeed, uniform convergence of ( u n ) and dominated 
convergence yield i 2 -convergence of (N Vj (u n )Vu). Hence, 

N(u n ) —>■ N(u), as n —> oo. 

By completeness of H 1 this implies N(u) G H 1 ( R + ) and also shows the continuity 
of N for the case to = 1. 

For the induction step from m to to +1 we may assume that the claim holds true 
for to > 1 and that Assumption A.2 holds for (m +1). Clearly, the assumption also 
holds for to and so N maps H m+1 continuously into H m by induction hypothesis. 
It thus remains to show that also o N maps H m+1 into H m . We decompose 
o N as in (A.l) and note that Assumption A.2 for to is also satisfied by 
and by 

(A.2) (x, y, z) y 3 (x, y, z) := -y(x, y)z. 
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Hence, by induction hypothesis the operators N x and Nj are continuous from 
H m ( R + ) d , resp. H m (M. + ) d+1 , into ), where Nj is the Nemytskii opera¬ 
tor defined by jlj, j = 1 Since also u M- Vu is continuous from H m+l into 

H m and Lemma A.l shows continuity of multiplication, (A.l) yields continuity of 
o N from 171+1 into H m , as claimed. □ 


Theorem A.7. Let p satisfy Assumptions A.2 and A.3 for some positive integer m. 
Then, N is Lipschitz continuous from bounded subsets of (H m (S. + )) d into 

Proof. We proceed as above, by induction on m G N. First, let m = 1 and u, 
v G H ^R+jR 7 *) with ||u|| ff i, |M| h i < r. By continuity of N we can assume, 
w. 1. o.g. u, v G iL 1 (R + ;R d ) D C°°(R>o;R d )- By Sobolev embeddings there exists 
a constant c s. t. ||u|| ,||t;|| < cr. By Assumption A.3, 

II N(u) — N(v )|| L 2 < L cr ||u — u||i 2 ( R+;R ij), 
ll-^x(w) — Aj,(n )|| L 2 < L cr ||u — I'||^,2 (m + ; K d ) > 

and for j = 1,..., d , 


(A.3) || N Vj (u)Vuj - N Vj (v)Vv 


Vj yu,j v LVJ j. vy. yuj v uj \\ L2 

< ||Ar % (u) 


Vuj - Vvj\\ L2 + \\N Vj (u) - N yj {v)\ 


IIL 2 


— Kj,cr || Vl/j S/Vj || L 2 + L cr \\llj Vj || L 2 , 


for Kj iCr 


ikj^y) 


< oo. Chain rule (A.l) then yields the 


®^P|y|<cr ®^Pa;£R 
assertion for m = 1. 

For the induction step we may assume that the theorem holds for fixed m and 
that Assumptions A.2 and A.3 are satisfied for m + 1. By induction hypothesis, 
N is Lipschitz on bounded sets from ff m (R + ) d into iL m (R + ) and thus, also from 
U m+1 (R + ) d into H m ( R + ). Hence, it suffices to show that ^ o A is Lipschitz on 
bounded sets from L/™(R_|_) d into ff m (R + ). To this end note that satisfies 
Assumptions A.2 and A.3 as well as pj , j = 1,..., d, defined in (A.2). By induction 
hypothesis, the operators N x and Nj, j = 1 defined in the proof of The¬ 

orem A.6, are Lipschitz on bounded sets. Again, approximation by elements in 
H m fl C m and (A.l) then show that the same holds true for ^oJV. □ 


Appendix B. The noise operator 

In this section we will study the operator-valued map C, defined in (2.13) by 

/ <t + (x,ui(x))(T c u;)(u 3+x) \ 

(C(u)w)(x) = I ct_(-x,u 2 (x))(T c w;)(w3 - x) I 

for u G 2?(A), w G U and x G R. We can reduce the problem to the operator 

(B.l) (u,x*) ct(.,r(.))T c (. + x*) 

for cr satisfying Assumption 2.3 and ( as in Assumption 2.6. Define the Nemytskii 
operator 

N a : H 2 (R + ) -a ff 2 (R+), u^ct(.,u(.)), 
which is Lipschitz on bounded sets by Theorem A.7. 

Lemma B.l. Multiplication is bilinear continuous from iL 2 (R + ) x Bf7C 2 (M>o) 
into H 2 ( R+). 
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Proof. By density of C n (~l H n in H n , neN, one can check that Leibniz formula 
holds for multiplication on H n x BUC n , so that 

D k {uf) = Y J ( k )^uf^ 

which is clearly square integrable for u £ H n ( R + ), / £ BUC n (WL> q) and k < n. In 
particular, for n = 2, 

k 

\D\uf){x)\ < K ||/|| C2 Y, &<*), k < 2 

3=0 

for some constant K, so that for some K. 

(B.2) KII^ (r+) <^II/II^IHIh- □ 

Lemma B.2. maps U into BUC 2 (WL). Moreover, T^w and its first two deriva¬ 
tives are Lipschitz continuous for all w £ U. 


Proof. First note that for all x, y £ R it holds that 


sup 

-l<e<l 


J^C W 0 + e,y) < j C {l+1) (x + e',y) de'+ ( M {x,y) 


i = 0 ,..., 2 


so that Lebesgue’s differentiation lemma gives 

(B.3) -^T^i)w(x) = T^ i+ i)w(x), x£R, w£U. 


Hence, it suffices to show that T^w £ BUC{ R) provided that (2.6) holds for i = 0 
and i = 1. For fixed w £ U and X\, x 2 £ R we directly get 


(B.4) 


iT^wt^Xx) - T c w(x 2 )\ < f [ \C(x 2 + e(xi - x 2 ),y)\\w(y)\ de dy \xi - x 2 | 

JR JO 


< sup||C , (a:,.)|| i 2 (R) INK |*i -x 2 \ 
zeR 


Here, we used fundamental theorem of calculus, Tonelli’s theorem and the Cauchy- 
Schwartz inequality. Hence, Tqw is globally Lipschitz and particularly uniformly 
continuous. Analogously, Cauchy-Schwartz inequality yields 


(B.5) sup|T f w;(*)| < sup ||C(z, .)|| L 2 (R) IMIu < oo. 


□ 


Remark B.3. From (B.5) we immediatly get Tf £ L(U, BUC 2 (R)). However, in 
general itself is not Hilbert-Schmidt. To get the Hilbert-Schmidt property we 
need the multiplication with N a as we will show in the next lemma. 

Lemma B.4. Let A be the Dirichlet Laplacian on L 2 (R + ). For u £ £>(A) and 
** £ R it holds that 

2 

IK K, ° -^C (’))) ll H S(U;X>(A)) <KMIa SU pX! |c w K0 

xeR i -0 

Remark B.5. This result immediately extends to C(u), because Assumption 2.3 and 
Lemma B.l assure N a (u) £ P(A) for all u £ V( A). Moreover, note that 

0 X °T{ = T Cx , 

where := f(x + .,.) satisfies Assumption 2.6, too. 
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Proof. Linearity and continuity in w follow directly from the construction and Re¬ 
mark B.3 and we are now interested in the Hilbert Schmidt norm. Without loss of 
generality, we can choose = 0. So denote by (e*,) an arbitrary CONS of U, then 

OO 

(B.6) |K(T C (.))|| 

HS(U,X>(A)) — Ell“'< I i e *)llh + ll A (“'( T c e ‘))llh. 
k =1 

and the first sum equals 

oo „ 

( B - 7 ) E / u ( x ) 2 (C( x i •)> e fe)i 2 ( R) 

fc=l ^ R + 

= [ u(x) 2 ||C(x, .)IIl 2 (r) < |I u IIi, 2 (r ) sup ||C■)IIl 2 ( m) . 

where we used Tonelli’s theorem and Parseval’s identity for the first equality. To 
bound the second sum we proceed on exactly the same way but first apply Leibnitz 
rule to get the second (weak) derivative 


(B.8) 


k —1 i=0 


EEL)/ £* u ( x ) <c (2 i \x,.),e k ) 2 L2( 


^E 

2=0 


dx l 


sup 

L 2 (R+) x gR 




2 

L 2 (R) 


2 

< 4 IMlL*(R+) sup |c w (a:, ■) 

„-_n 


2 

L 2 (R) ‘ 


By equivalence of ||.|| ff2 and ||.[| A the result follows. 


□ 


To show the main result of this appendix, we just need to combine the previous 
lemmas. 


Theorem B.6. The map C : 'D(A) —> HS(f7, D(A)) is Lipschitz continuous on 
bounded sets. 


Proof. By the structure of A and £ 2 , it suffices to show the property for the operator 
defined in (B.l). Assumption 2.3 yields N a (T>(A)) C 2?(A) and we can apply 
Lemma B.4. For u, u £ 2?(A), x*, y* £ M and writing Q t z(x, y) := f(z + x,y) — 
f(z + x), it holds that 

\\N a (u)(9 x T c ) - NriuWyT^ 

< m*(u) ~ N a (u)) ■ (0*r c .)|| HS + 11^(5) • (0 o T^ y .)\\ m 

2 

< R c WN^u) - N a (u)\\ A + K ||7V CT (n)|| A sup^ 

zGR • n 



A computation similar to (B.4) shows 


sup 

zGR 




y I sup 


zei 


C (i+1) (*,0 


2 

L 2 ’ 


Finally, we put everything together and use that on bounded sets N a is Lipschitz, 
and thus bounded, to get the assertion. □ 
























A STEFAN-TYPE STOCHASTIC MOVING BOUNDARY PROBLEM 


33 


[AZ90] 

[BDP02] 

[BHQ14] 

[BMS05] 

[CKS14] 

[DBMB14] 

[DPZ14] 

[DZ92] 

[EN001 

[EPS03] 

[FP79] 

[Gri66] 

[Hen81] 

[Kal02] 

[KMS10] 

[KS12] 

[Kur66] 

[KZS12] 

[LL07] 

[LM72] 

[LPS13] 


References 

J. Appell and P. P. Zabrejko. Nonlinear superposition operators , vol¬ 
ume 95 of Cambridge Tracts in Mathematics. Cambridge University 
Press, Cambridge, 1990. 

V. Barbu and G. Da Prato. The two phase stochastic Stefan problem. 
Probability Theory and Related Fields , 124(4):544-560, 2002. 

C. Bayer, U. Horst, and J. Qiu. A Functional Limit Theorem for 
Limit Order Books with State Dependent Price Dynamics. ArXiv e- 
print:1405.5230, May 2014. 

Z. Brzezniak, B. Maslowski, and J. Seidler. Stochastic nonlinear beam 
equations. Probability Theory and Related Fields , 132(1):119 149, 2005. 
R. Cont, A. Kukanov, and S. Stoikov. The price impact of order book 
events. Journal of Financial Econometrics, 12(l):47-88, 2014. 

J. Donier, J. Bonart, I. Mastromatteo, and J.-P. Bouchaud. A fully 
consistent, minimal model for non-linear market impact. ArXiv e- 
print: 1412.0141 , November 2014. 

G. Da Prato and J. Zabczyk. Stochastic equations in infinite dimen¬ 
sions, volume 152 of Encyclopedia of Mathematics and its Applications. 
Cambridge University Press, Cambridge, second edition, 2014. 

G. Da Prato and J. Zabczyk. A note on stochastic convolution. Sto¬ 
chastic Analysis and Applications, 10(2):143-153, 1992. 

K. -J. Engel and R. Nagel. One-parameter semigroups for linear evolu¬ 
tion equations, volume 194 of Graduate Texts in Mathematics. Springer- 
Verlag, New York, 2000. 

J. Escher, J. Priiss, and G. Simonett. Analytic solutions for a Ste¬ 
fan problem with Gibbs-Thomson correction. J. Reine Angew. Math., 
563:1-52, 2003. 

A. Fasano and M. Primicerio. Free boundary problems for nonlinear 
parabolic equations with nonlinear free boundary conditions. J. Math. 
Anal. Appl, 72(l):247-273, 1979. 

P. Grisvard. Commutativite de deux foncteurs d’interpolation et appli¬ 
cations. J. Math. Pures Appl. (9), 45:207-290, 1966. 

D. Henry. Geometric Theory of Semilinear Parabolic Equations, volume 
840 of Lecture Notes in Mathematics. Springer Berlin Heidelberg, 1981. 
O. Kallenberg. Foundations of Modern Probability. Applied probability. 
Springer, 2002. 

K. Kim, C. Mueller, and R. B. Sowers. A stochastic moving boundary 
value problem. Illinois J. Math., 54(3):927 962, 2010. 

K. Kim and R. B. Sowers. Numerical analysis of the stochastic moving 
boundary problem. Stochastic Analysis and Applications, 30(6):963- 
996, 2012. 

K. Kuratowski. Topology, volume 1. Acad. Pr., 1966. 

Ku. Kim, Z. Zheng, and R. B. Sowers. A Stochastic Stefan Problem. 
Journal of Theoretical Probability , 25:1040-1080, 2012. 

J.-M. Lasry and P.-L. Lions. Mean field games. Jpn. J. Math., 2(1):229- 
260, 2007. 

J.-L. Lions and E. Magenes. Non-homogeneous boundary value problems 
and applications - 1. Springer, Springer, 1972. 

A. Lipton, U. Pesavento, and M. G Sotiropoulos. Trade arrival dynam¬ 
ics and quote imbalance in a limit order book. ArXiv e-print: 1312.0514, 
December 2013. 


34 MARTIN KELLER-RESSEL AND MARVIN S. MULLER 

[Lun95] A. Lunardi. Analytic Semigroups and Optimal Regularity in Parabolic 
Problems. Progress in Nonlinear Differential Equations and Their Ap¬ 
plications. Birkhauser Basel, 1995. 

[Lun04] A. Lunardi. An Introduction to Parabolic Moving Boundary Problems. 

In Mimmo Iannelli, Rainer Nagel, and Susanna Piazzera, editors, Func¬ 
tional Analytic Methods for Evolution Equations , volume 1855 of Lec¬ 
ture Notes in Mathematics, pages 371-399. Springer Berlin Heidelberg, 
2004. 

[Lun09] A. Lunardi. Interpolation theory. Appunti. Scuola Normale Superiore 
di Pisa (Nuova Serie). [Lecture Notes. Scuola Normale Superiore di Pisa 
(New Series)]. Edizioni della Normale, Pisa, second edition, 2009. 

[MTB14] I. Mastromatteo, B. Toth, and J.-P. Bouchaud. Anomalous impact in 
reaction-diffusion financial models. Phys. Rev. Lett., 113:268701, Dec 
2014. 

[Paz92] A. Pazy. Semigroups of Linear Operators and Applications to Partial 
Differential Equations. Number 44 in Applied Mathematical Sciences. 
Springer, 1992. 

[PR07] C. Prevot and M. Rockner. A Concise Course on Stochastic Partial 
Differential Equations. Springer, 2007. 

[PSS07] J. Priiss, J. Saal, and G. Simonett. Existence of analytic solutions for 
the classical Stefan problem. Math. Ann., 338(3):703-755, 2007. 

[Ste88] J. Stefan. Uber die Theorie der Eisbildung, insbesondere iiber die Eis- 
bildung im Polarmeere. Wien. Ber. XCVIII, Abt. 2a (965-983), 1888. 

[Val85] T. Valent. A property of multiplication in Sobolev spaces. Some appli¬ 
cations. Rend. Sem. Mat. Univ. Padova, 74:63-73, 1985. 

[Val88] T. Valent. Boundary Value Problems of Finite Elasticity: Local The¬ 
orems on Existence, Uniqueness, and Analytic Dependence on Data, 
volume 31 of Springer Tracts in Natural Philosophy. Springer New 
York, 1988. 

[vNVW12] J. van Neerven, M. Veraar, and L. Weis. Maximal L p -regularity for 
stochastic evolution equations. SIAM J. Math. Anal., 44(3):1372-1414, 
2012. 

[Vui93] C. Vuik. Some historical notes about the stefan problem, 1993. 

[Wer07] D. Werner. Funktionalanalysis. Springer-Lehrbuch. Springer Berlin 
Heidelberg, 2007. 

[Zhel2] Z. Zheng. Stochastic Stefan Problems: Existence, Uniqueness and Mod¬ 
eling of Market Limit Orders. PhD thesis, Graduate College of the 
University of Illinois at Urbana-Champaign, 2012. 

Institut fur Math. Stochastik, TU Dresden, Germany 

E-mail address, M. Keller-Ressel: martin.keller-ressel8tu-dresdeii.de 


E-mail address, M. S. Muller: marvin.s.mueller@tu-dresden.de 


